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Abstract 

The quasi-Gaudin algebra was introduced to construct integrable systems which are only quasi- 
exactly solvable. Using a suitable representation of the quasi-Gaudin algebra, we obtain a class of 
bosonic models which exhibit this curious property. These models have the notable feature that 
they do not preserve (7(1) symmetry, which is typically associated to a non-conservation of particle 
number. An exact solution for the eigenvalues within the quasi-exactly solvable sector is obtained 
via the algebraic Bethe ansatz formalism. 

PACS Numbers: 02.30.Ik, 03.65. Fd, 05.30.Jp. 

1 Introduction 

In [Tl[2] Ushveridze proposed a method for studying quasi-exactly solvable (QES) systems ^^S, from the 
prespective of integrable systems and the Quantum Inverse Scattering Method (QISM) [6]. The approach, 
which is called the partial algebraic Bethe ansatz (ABA) , relies on deforming the Yang-Baxter algebra in 
such a way that it retains most of the features required for the QISM but leads to generating functions of 
integrable systems which are only QES. This deformation of the Yang-Baxter algebras led to new classes 
of hitherto unknown algebras. A limiting case is the (rational) quasi-Gaudin algebra which will be the 
focus of this study. 

Exactly solvable models have found many successes in various branches of physics and mathematics. 
Over recent years they have continued to find new applications in diverse fields such as Bose-Einstein 
condensates and degenerate Fermi gases, quantum optics, superconductivity, and nuclear pairing among 
other things e.g. [7Hl4]. There has also been significant interest in QES models, with new applications 
of these being found in problems relating to matrix product states [15], and in dissipative systems [16]. 
However by comparison the partial ABA approach seems to have received little attention and remains 
essentially undeveloped. An appealing property of the partial ABA is that it provides us with a con- 
structive algebraic approach for obtaining QES models which have multiple degrees of freedom. 



One particular aspect of the ABA which bears some relevance to our present exposition is the study 
of quantum integrable models which do not preserve U{1) symmetry. Such models are interesting for a 
number of reasons. In the context of spin-boson Hamiltonians of the Tavis-Cummings form, these models 
correspond to physical systems without the rotating wave approximation. Diagonalisation of such models 
is a somewhat complicated affair within the ABA method due to the lack of reference states, often re- 
quiring the use of functional Bethe ansatz or Sklyanin's separation of variable technique [51ITT|. Non U{1) 
preserving models are also relevant to the study of open quantum systems whereby the U{1) symmetry 
is broken due to coupling to an environment. An example of this is found in the spin-boson Hamilto- 
nian of Leggett et al [17] which has found applications ranging from quantum-state engineering [18] to 
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biomolecular systems |19j . 




In the present paper, we will study QES bosonic models descending from suitable realisations of the 
quasi-Gaudin algebra. It will be shown that such models correspond to an extension of the su(l, 1) Dicke 
Hamiltonian |20) by the addition of U(l) symmetry-breaking terms. The Hamiltonian can be written as 

H = Ho + Hi (1.1) 

with Hq and Hi refering to the Dicke Hamiltonian and the U(l) symmetry-breaking component respec- 
tively. Explicitly, they have the form 

Ho - wNb 

Hi = g(^{b + b^)\^n + r-Y,S!j^b^b'-ib^)Hj. (1.2) 

Here Ni,,b,b'' are standard bosonic operators, f is a representation dependent parameter, wo,ei,g are 
free parameters, n is an integer and Sf'^ are either single-mode or double-mode representations of su{l, 1) 
generators (refer to equations (|3.22p and (13.23^ below). The Hamiltonian Hi may be interpreted as a 
coupling of the su(l, 1) Dicke model to an external system. 

Our paper is structured as follows. In Section 2 we will briefly review the partial ABA method of 
obtaining quasi-exact solutions for models associated to the quasi-Gaudin algebra. In Section 3 we will 
use a suitable representation of the quasi-Gaudin algebra to obtain the integrable bosonic model (1.1). We 
then derive the Partial ABA solution of the Hamiltonian and discuss aspects relating to the quasi-exact 
solvability. Finally in Section 4 we summarise our results and discuss possible future lines of work. 

2 Quasi-Gaudin Algebra and Bethe Ansatz Solution 

Let us first introduce the rational (rank 1) Gaudin algebra and the associated abstract, integrable models 
before defining its quasi counterpart. The rational Gaudin model is a parameter-dependent infinite- 
dimensional Lie algebra satisfying the following commutation relations: 

S'{X)S^fi)~ S'{fi)S'iX) = 0, 
S''{X)S^{fi)~S^{fi)S%X) = 



S-{X)S+{fi) - S+{n)S-'{X) = 2 



/i — A 
5-(A)-5-(m) 



/i — A 

whereby A and /i are complex spectral parameters. From these relations, it can be shown that 

H{X) = S'^(A)S'^(A) - is'+(A)5-(A) - ^S-{X)S+{X) (2.3) 

satisfies the following commutation relations 

[H(A),iJ(M)]=0 (2.4) 

and therefore acts as a generator of commuting operators in an abstract integrable system. Assuming the 
existence of a suitable reference state, the spectrum of H{X) can be obtained via the standard ABA [2]. 
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Analogous to the Gaudin algebra is the so-called quasi-Gaudin algebra. It is defined by the following 
parameter dependent set of relations [TJ[5] 

s^Ws^iii) - s^Ms^^W = 0, 

5±±i(A)5±(m)-^^±i(a*)5±(A) = 0, 

jl — A 

S-+,{X)S+{i,)-S+_MSnW = 2 ^"^^^~f (2.5) 

/i A 



whereby n is an integer and A, /i are complex parameters. While \2.h\ appears to be similar to the 
Gaudin algebra, we stress that there are important qualitative difference between the two. Importantly, 
note that (|2.5p do not define commutation relations and are therefore not Lie algebraic relations. Despite 
looking somewhat arbitary, the quasi-Gaudin algebra can be understood as a grading deformation on the 
original Gaudin algebra. We refer the reader to [2^ for a more detailed discussion. 

Similar to the Gaudin algebra, there exists a generating function of commuting operators for the 
quasi-Gaudin algebra. It has the form 

i/„(A) = S:X>)S:,{\) - i^-+i(A)5+(A) - \s+_,{X)S-{^) (2.6) 

and can be shown to form a commutative family with respect to the spectral parameters, i.e. 

[i/„(A),i/„(Ai)] =0. (2.7) 

Note that the commutation relation (|2.7I) does not extend to the general case where i?„(A) and H,n{^) 
have different integer values of n and m. This is due to the lack of a defining relations between elements 
of the algebra with arbitrary integer indexes. The ABA solution for the generating function Hn(X) of 
the quasi-Gaudin algebra has been obtained in [ll[2]. As wtih the standard Gaudin algebra, the ABA 
diagonalisation of i?„(A) works if the representation of (|2.5I) supports a reference state |0), viz. 

5o^(A)|0)=/(A)|0) , %(A)|0)-0 (2.8) 

The Bethe vector is given by 

■■■fin)^ 5+_i(/.„)5+_2(Mn-l) • ■■S+{fli)\0). (2.9) 

By successively applying the following relation 

i?„(A)5^_i(Ai„) = S'„_i(Ai„)i?„-i(A) + 2 (2.10) 

we can shift the operator i/„(A) towards the right of the product of S^{iJ,i+i) operators on the right- 
hand side of (j2.9p . so that we finally have i/„(A) acting on the reference state. After having completed 
this procedure, we perform the same operation for the various S'?(A), S"? (/ii+i) that were generated as a 
byproduct of shifting the i?„(A) through the product of the S^{iJ.i+i). The final form is given by 

H{\)tp(pi,- --fin) = A{X)'ll;{fii, ■■■Hn) +2^B{fli)^{ni,- ■ ■ X,fit+1, ■ ■ ■ , fJ.n) (2.11) 

i 

whereby 

AiX) = fiXr + f'W + 2±J^ + 2±^j:-^^ 
B{^^^) = f{fi^) + Y,^—- (2-12) 
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By requiring that the unwanted terms vanish we obtain the fohowing Bethe ansatz equations: 

n ^ 

-^— + f{f^i)^0, z = l,2,...,n (2.13) 



En{\) = f{\) + /'(A) + 2 f^^^ f^^'^\ (2.14) 



with the eigenvalue for -ff„(A) given by 

As a proof of existence, an explict representation for (|2.5p is provided in [Hll]: 

S-{X) ^ 5-(A)+'^'"^'+" 



A-c 



A-c 

5+(A) = 5+(A) + ^l-^±^^. (2.15) 

A — c 

with c and c? as free parameters, S^'^(X) are generators of the Gaudin algebra, and and are defined 
as 

S-^ = lim AS'^(A), S-^IO) = /"|0). (2.16) 
In terms of this realisation, the generating function Hn{X) takes the form 

i7„(A) = S'{X)S''{\)-^S-{X)S+{X)-^S+{X)S-{X) 

2S''{X){n + d+f'-S')-S-{X){n + 2d+f'-S')-S+{X)(n + f'-S') 



A-c 

1 



4(A-c)^- (2-^^) 

It can be seen that the condition of hermiticity for (|2.17p is satisfied when d = 1/2 and the representation 
for the Gaudin algebra is unitary, i.e. satisfying the condition 

S+{X)^ = 5-(A), 5^(A)t = S^iX). (2.18) 

3 Bosonic Representations of the Quasi-Gaudin Algebra 

The quasi-Gaudin algebra of the form (|2.15l) admits mixed representations, consisting of su(l, 1) algebras 
and the Heisenberg algebra, with the following form: 

5-(A) = ^ + ' r-N.-E.S! + n 



9 A - ej A-c 



i=l ' ' ^ 



^-^ A — A-c 



a ^-^ A — € i X — c 

i—l 

(3.19) 

The S^'^ and {Ni,,b,b^} are respectively the sm(1, 1) and Heisenberg algebras, which obey the commu- 
tation relations 

[Si , S^] — ±SfSij , [S^ , 5^] —23^ 6ij 

[m, &t] ^ 5t ^ [N^^ b]^~b, [b, 6t] = 1 (3.20) 
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and and are defined as 

m 

= Y,S- +Ni„ 5"|0) = /^|0). (3.21) 

i=l 

We note here that our definition for differs from that of (|2.16p as the prior definition is divergent for 
this particular reahsation. 

The su(l,l) algebras has two bosonic operator reahsations. The first is given by the single-mode 
representation, whereby 

-,t„. 1 AT 1 lJ\'i „2 



ojo, 1 ^ 1 + ^ K)^ 

2 4 2 4'' 2 



St = ^ + 1 = ^ + 1. S+^^ , Sr^^. (3.22) 



The second one is given by the two-mode representation, 



St = ^[aS + clc,)+- = ^^^^^-^ + -, S: = alcl, 5" = a,c,. (3.23) 

There are multiple reference states for both bosonic realisations. For the single-mode realisation, there 
are finitely many of them. We can express them as 

m 

|0,{/}) = n(«l)''|0) ' ^. = Oorl (3.24) 

1=1 

where {1} is a shorthand notation for the set {/i, • • • Im} and 

5^lO,W) = nO,W)= (^f;| + i^ |0,W). (3.25) 

For the two-mode reahsation, there are infinitely many reference states. Without loss of generality we 
can write them as 

m 

|0,{O) = nK^)''|0)' ^. = 0,l,2,--- (3.26) 

1=1 

with 

^10, {I}) ^ no, {I}) ^ I + 1^ |0, {I}). (3.27) 

It can be seen that each reference state corresponds to a distinct eigcnfunction of the Casimir operators 
for the sm(1, 1) generators S^'^ . As the sm(1, 1) Casimir operators acts as central elements with respect 
to (|3.22p and (|3.23p . we can use Schur's lemma to deduce that each reference state gives rise to a distinct 
irreducible representation. 



4 Quasi-Exactly Solvable Hamiltonians 

We now consider the generating function Hn{\) of the quasi-Gaudin algebra obtained from the represen- 
tation (|3.19p . Assuming ^ ej, it can be seen that 




(4.28) 
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with 



m ^ 



+- 



c 



2 ^ (n+l + .r-E"i>gf-A^b ) - + r - S^) - ^Sr{n + 1 + - S^) 

tic ^ g 2^ 



=1 (c-Ci) 



= - i (5r^+ + 5+5r) • (4.29) 

From (|4.28p and the commutation relation (|2.7p . it follows that Hi^c,Ki.c form a set of mutually com- 
muting operators. By considering the following linear combination H ^ T + Hi + He and setting the 
coefficient c = 0, we obtained the desired bosonic hamiltonian. For the single-mode representations, we 
have 



m m 



{n + n{b^ + 6) - (6 + 6t) ^ ^ _ 5t52 _ (^t)2^ j (4 30) 



where 



T = w n , 

2 — ^2 



j=i 



For the two-mode representations, we obtain 

m m 
2=1 i=l 

+5 (^(" + r + 6) - (6 + 5t) ^ ^ - - (6t)26^ (4.31) 

with 

We note that for the case when m = 1, the models correspond to quasi-exactly solvable extensions for 
atom-molecule BEC models contained |21| . 

The eigenvalues for the Hamiltonians can be extracted from the Bethe ansatz solution of (|4.28p : 

E,,{\) = /2(A) + /'(A)+2f]M^-/i^. (4.32) 

This is done by evaluating the residues of the poles and c. Doing so yields 

whereby sf = {2li + l)/4 for the single-mode representations and sf = {k + l)/2 for the two-mode 
representations . 
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We now examine the quasi-exactly solvable nature of the Hamiltonians in more detail. For the sake 
of clarity, we shall only consider the Hamiltonian with the single- mode bosonic representation (j4.30p . as 
results for the two-mode representation will follow analogously. It is straightfoward to see that (|4.30|) 
acts on an infinite-dimensional Hilbert space V span by the following basis states 

V = span{(6t)'o („t)ii . . . jg)} ^ span{|/o, • • • , /„.)}, h £ Z+ . (4.34) 

In order to identify the invariant subspace which characterises the quasi-exact solvability of the Hamil- 
tonian, let us write the Hamiltonian as 

Hg = Ho + H^+H+ (4.35) 
whereby we have introduced a grading structure on the Hamiltonian through setting 

rn / rn \ 

Ho = «;7V,+^e,7V,. +5K]6(al)2-h6ta2 + (^ + j~.)(^t+5)^ 
j=i \i=i / 



i=l 

2 



= (n + /^)6-6^^ -6+62 (4.36) 
1=1 

The assigned grading of ±,0 is determined by the commutation relations of i?±,o with the C/(l) charge 
5^ = A^b + Elli(2A^a, +l)/4: 

[S'^i^°]=0 , [S\H+]^H+ , [S''-,H-]^ -H-. (4.37) 

In light of these relations, we may decompose V into a direct sum of eigenspace Vi^p of the U{1) charge 
and the Casimir operators of the sm(1, 1) algebra Ki = S^{Sf — 1) — S~ , i.e. 

^=0^^M- (4-38) 
Explicitly, the subspace can be written as 

y,,{p|=span{(6t)'o(4)2'i+fi...(a1„)2'".+f™|0)}, = or 1. (4.39) 

3=0 

It can also be verified that 



From the commutation relations (j4.37p . we therefore have 

H+Vi^{p] C Vi+i^{p-}, HoVi^{p] C Vi^{p}, H^Vi^{pT^ C yj„i_|p}. (4.41) 

The QES property of the Hamiltonian arises from the fact that for given integer value of n and — 
Y,, ik + 1) /4, we have 

i?+F„,{,} = {0}. 

As a result, the Hamiltonian leaves the following subspace invariant: 

n 

^QES^0^.,m- (4-42) 

1=0 

We can indeed verify that the Bethe vectors lie within this invariant subspace, by expanding the eigen- 
vectors (12.91) explicitly. It would be interesting to examine the possibility of obtaining exact solutions 
outside of this sector. 
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5 Conclusion 



We've investigated a class of QES, integrable multi-mode bosonic models using the quasi-Gaudin alge- 
bra. We see that such models are obtained via a mixed representation consisting of commuting copies 

of sii(l,l), and the Heisenberg algebra. Integrable Hamiltonians were extracted from the generating 
function of commuting operators. A notable feature was that the QES Hamiltonians we obtain do not 
preserve U{1) symmetry. We identified the QES sector of the Hamiltonian as the direct sum of the 
eigensubspaces of the U{1) charge with eigenvalues no greater than n. 

The ABA method leads to partial solutions of the Hamiltonians we've considered. Given the inte- 
grability of the Hamiltonian, in the sense that i?„(A) acts as a generator of conserved operators, it would 
be interesting to explore the possibility of obtaining the entire spectrum via some other techniques. The 
dominating experience is that integrability and exact solvabilty go hand-in-hand. It is not apparent for 
these Hamiltonians whether the full spectrum is potentially accessible. 

Finally wc note that due to the constraint arising from imposing hermiticity on the generating function 
Hn{\), the quasi-Gaudin formalism is at present limited to cases based on underlying unitary represen- 
tations of su(l, 1), or the Heisenberg algebra. It would be of interest to obtain representations of the 
quasi-Gaudin algebra based on non-unitary (in particular finite-dimensional) representations of su{l, 1), 
which are also able to accomodate hermitian Hamiltonians. 
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